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We investigate the microwave absorption of a pinned, two-dimensional Wigner crystal in a strong 
magnetic field at finite temperatures. Using a model of a uniform commensurate pinning potential, 
we analyze thermal broadening of the electromagnetic absorption resonance. Surprisingly, we find 
that the pinning resonance peak should remain sharp even when the temperature is comparable or 
greater than the peak frequency. This result agrees qualitatively with recent experimental observa- 
tions of the ac conductivity in two-dimensional hole systems in a magnetically induced insulating 
state. It is shown, in analogy with Kohn's theorem, that the electron-electron interaction does not 
affect the response of a harmonically pinned Wigner crystal to a spatially uniform external field at 
any temperature. We thus focus on anharmonicity in the pinning potential as a source of broad- 
ening. Using a expansion technique, we show that the broadening is introduced through the 
self-energy corrections to the magnetophonon Green's functions. 

PACS numbers: 73.40.Hm, 73.20. Dx 



I. INTRODUCTION 



It is now well appreciated that in certain situations, the translational symmetry of an electron gas becomes broken 
at low temperatures giving rise to an electron solid state called the Wigner crystal (WC).El Among other conditions, 
low density and low disorder are crucial for this novel state to be realized; the inter-electron distance has to be large 
enough for the Coulomb interaction to dominate the quantum zero-point fluctuations, and the disorder must be weak 
enough not to destroy the long- (or quasi-long-) range order. One of the systems that satisfy these conditions best 
is the two-dimensional (2D) electron gas in the GaAs/AlGaAs heterostructure, in which density and disorder can 
be controlled relatively easily. Furthermore, if a strong magnetic field is applied perpendicular to the 2D plaae, the 
cyclotron motion quenches zero-point fluctuations and the low-density condition is even more easily achieved.B 

Although direct observations of the positional order from scattering measurements, for example, have been unattain- 
able, there have been considerable, efforts to prove the existeroce of the WC through more or less indirect methods 
such as transport measurementaJcl and photoluminescenceOiij One of the most convincing evidences so far is that 
the 2D systems become insulating when the fllling factor ly — nhc/eB is low: ^1/5 for electron systemaJ and ^1/3 
for hole systems.l3 This insulating behavior is generally accepted as a result of pinning of the WC due to impurities. 
The impurity potential is also supposed to provide a restoring force in the ac electromagnetic response of the WC. 
If the frequency of the external driving force matches the natural frequency of the pinning mode, a resonance should 
occur. rj |-| 

Recent experiments on 2D hole systemsQij have revealed a resonance structure in the ac absorption spectrum in 
the low fllling factor regime. However, some features of the spectrum were qualitatively different from the predictions 
of the previous calculations based on charge-density-wave (CDW) models;ll3E3 namely, (i) the peak frequency ujp 
increased with the magnetic fleld B, and (ii) the quality factor Q — ujp/Auj was found to be as large as 30 in some 
data.cl The CDW calculations had predicted that ujp would decrease with B and the peak would be much broader 

These puzzling experimental flndings were later examined by the authors of Refs. Oand nM Based on weak disorder 
models at T = 0, they found that (i) was explained if the disorder potential was allowed to vary within a magnetic 
length Iq — y^hc/eB. This short-length-scale physics is missing in the previous CDW calculations since they assume 
that the disorder potential can vary only in a length scale longer than the lattice periodicity a. (In the experiments 
in Refs. 0, ^, and ^, a > Iq.) More specifically, in Ref. |8[ it was suggested that the roughness of the GaAs/AlGaAs 
interface could serve as a source of such a short-length-scale disorder potential. A "pit" on the interface, which is 
typically several atoms wide, can trap an electron by allowing it to stay closer to the positively charged donor plane 
and gain electrostatic energy. In a strong magnetic field, one can project the many-body wave function onto the 
lowest Landau level (LLL). Assuming a ^ Iq, one may ignore exchange energy and describe the many-body ground 
state as a Hartree-type product, of single electron wave packets of size ~ Iq, which are, in the absence of disorder, 
located at crystal lattice sites.E^I Suppose an electron wave packet is situated right on a pit. Since its size decreases 
with increasing B, the probability for the electron to stay inside the pit subsequently increases, provided the pit size 
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is smaller than Iq- Since LUp is essentially given by the expectation value of the pinning potential strength, the above 
cartoon picture explains how Up grows with B. In addition, as described below, the long-range Coulomb interaction 
suppresses low-energy collective excitations and leads to a sharp resonance. 

At finite temperatures, however, thermal fluctuations should broaden the resonance peak. One important observa- 
tion is that the peaks measured in experiments keep getting narrower down to the lowest T measured. This suggests 
that the main source of broadening is thermal fluctuations. Interestingly, the experiments also find that Auj <C T 
even though T ^ ujp. Thus, how the thermal broadening comes about is a highly non-trivial question. We address 
the question in this paper. As the thermal broadening is the main issue, we will simplify the disorder model and 
assume that the pinning centers themselves form exactly the same lattice as the WC. We also assume that the im- 
purity potential has the same strength at all sites.E3 Several different disorder potentials — incjluding the uniform 
commensurate model we will use in this work — have been studied at T = in a previous study,E3 but it was found 
that the details of disorder do not alter the resonance structure qualitatively. We believe the thermal behavior is not 
sensitive to details of disorder, either. 

Our result can be summarized as follows. In order to understand low, but finite temperature behavior, one needs 
to investigate the low-energy excitations and their interactions. The longitudinal and the transverse phonon modes of 
the WCLmp. together in a strong magnetic field and two new normal modes, magnetophonons and magnetoplasmons, 
emergeEJo We focus on the magnetophonons because the energy of the other mode, magnetoplasmons, is essentially 
given by the cyclotron energy hujc = heB/mc, which is a few orders of magnitude greater than T and ujp. In a theory 
that treats displacement of electrons from lattice sites harmonically, magnetophonons are independent and produces 
only a delta-function peak in the ac absorption spectrum. The interaction between magnetophonons comes about 
from anharmonicity of the Hamiltonian. There are two sources of anharmonicity in our model: Coulomb interaction 
and pinning potential. rSince the external ac field has a much larger wave length (> 1 m) than the size of the system 
(typically ^ 10^^ m),Ll it is safe to treat it as a uniform field that couples only to the center-of-mass degree of 
freedom. If the pinning potential is completely harmonic, we find that the anharmonicity in the Coulomb interaction 
does not contribute to the absorption spectrum, because the interaction depends only on the relative position of 
electrons. This is analogous to the Kohn's theoremjHJ which states that the electron-electron interaction does not 
affect the cyclotron frequency in a disorder-free system. The theorem remains true even in a harmonically pinned 
WC, Therefore, anharmonicity of the pinning potential is crucial for understanding thermal broadening. 

As will be shown below, magnetophonon modes are very closely related to angular momentum excitations. It is 
therefore natural to exploit analogies to 2D spin lattice systems. In fact, written in magnetophonon creation and 
annihilation operators, our iwd^l Hamiltonian resembles that of an SU(N) Schwinger boson formulation of a 2D 
quantum (anti-)ferromagnet.tj E3 Further utilizing this analogy, we use a 1/7V expansion technique that has proved 
useful in spin systems. Not only does it provide a systematic way of collecting important processes in Feynman 
diagrams, but also its lowest-order mean-field (MF) solution already contains terms that are higher order in coupling 
constant in a standard diagrammatic many-body calculation. Importantly, the position of absorption resonance peak 
is mostly determined at the MF level. 

The Coulomb interaction and the pinning potential affect the spectrum of the magnetophonons significantly. With- 
out pinning, the_dispersion relation is well knownti tp take the form £q cx |q|'^/^ in the long wave length limit, calculated 
both classicallyEll and quantum mechanicallyHEj. In the uniform commensurate pinning model, it is modified to 
£q~ V + Ws|q|, where V represents the strength of the pinning potential and Vg is the sound velocity of the magne- 
tophonon mode. The new spectrum is not only gapped, but also makes the density of states (DOS) vanish linearly 
in the low-energy limit. It turns out that the slope of the low energy DOS is so small that there are extremely few 
low-energy excitation modes about cOp. As a consequence, we find that the thermal fluctuations are substantially 
suppressed and the peak remains sharp. 

This paper is organized as follows. The magnetophonon creation and annihilation operators are constructed out 
of displacement operators in the LLL in Sec. |l| We derive the ac conductivity in terms of magnetophonon Green's 
functions in Sec. [II. We then introduce in Sec. IV, the Hamiltonian of the WC, and derive in Sec. 0, a generalized 



Kohn's theorem that explains why the anharmonicity in the Coulomb interaction does not contribute to thermal 
shift or broadening of the absorption spectrum in a harmonic pinnin g mo del. Sec. ^ is devoted to 1/N expansion 
calc ulatio ns of the Green's functions. We present the results in Sec. VH, and finally conclude with a discussion in 
Sec. VIH. Details of some calculations are given in Appendices. 



II. MAGNETOPHONONS 



Before we write down the Hamiltonian in a second-quantized form, it is helpful to see how the Hilbert space is 
spanned by the eigenfunctions of magnetoplasmon and magnetophonon number operators. In the circular gauge, 



2 



A = (1/2)B X r, a 2D single electron wavefunction in the n-th Landau level with angular momentum quantum 
number m takes the form 



/ /o/ \2n+m _l^r 

ijnmir) = ^ e^9^"9^"+"e ^, (2.1) 



where z — x + iy, dz = {d/dx — id/dy)/2, dz* — {d/dx + id/dy)/2, and lo = ^Jhc/eB is the magnetic length. This 
wavefunction spans the entire single particle Hilbert space. It is now useful to define ladder operators that raise or 
lower n and to. We define the Landau level raising operator 
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and the angular momentum raising operator 



aVnm OC ■>pn+l,m-l, aipOm = 0, 

[a,a^] = [b,b^^l, 

[a, b] = [flt , b1] = [a, b1] = [at , 5] = 0. (2.4) 

In what follows, the above ladder operators will be defined at each WC lattice site. Then, a Fourier transformed 
operator will serve as a magnetoplasmon creation operator, and 6jj as a magnetophonon creation operator. We 
will use the latter terms to refer to the and b^ operators from this point on. 

In the low filling factor limit, the Landau gap htOc is large, so one can safely truncate the Hilbert space and work in 
the LLL only. Furthermore, the ground state of the WC in this limit is believed to be well represented by a product 
of a collection of zero angular momentum single particle wave functions, each of which is centered at a lattice site.E£l 
This is analogous to a Hartree wave function. Although the single particle wave functions are not orthogonal, their 
overlap is exponentially small if a- Thus, one may ignore exchange energy correction and the wave functions do 
not have to be antisymmetrized. In other words, electrons are distinguishable because each one of them is strongly 
localized at a lattice site far from the others. For higher angular momentum states, the size of wave packets grow with 
TO. However, the above argument should be still valid for the first few higher angular momentum states. Therefore, 
we assume that the Hilbert space is given by a direct product of single particle Hilbert spaces spanned by angular 
momentum states in the LLL. 

Due to the strong Coulomb repulsion, we assume that there is only one electron per site. This allows us to label 
each many-body energy eigenstate by a collection of Ng angular momentum quantum numbers or magnetophonon 
occupation numbers, where Ns is the total number of sites. As will be shown below, the above statement is not 
exactly true, because in reality, different angular momentum states mix together due to quantum fluctuations. This 
is demonstrated in our Hamiltonian as terms that do not conserve the magnetophonon number. Nonetheless, higher 
angular momentum mixing is small at low temperatures and the above argument should be a good qualitative 
description of the system. j_. 

The many-body wavefunction is written asE3 

^{{n})^llij„,,irf,n,), (2.5) 

i 

where the single particle wavefunction ipmiiju'Ri) is localized at a lattice site Rj and has an angular momentum 
quantum number m^. More specifically, it may be written as 

1 . ^ - 



z 



It is straightforward to show that 



where Zi = Rf + iR\- Note that the last exponential is a pure phase and is needed to keep the wavefunction in the 



LLL . We will assume that {Ri} forms a triangular lattice in the ground state. The ladder operators in Eqs. (2.2) and 



(2.3) are modified at each lattice site and may be written as 
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(2J) 



where Q = Zi — Zi is the complex notation of the displacement of the «th electron from its site. 
In terms of the above ladder operators, the original electron position is written as 

X, = i?f + -^(6, + 6l + a, +a|), (2.9) 
v2 

y^ = Ry + ^{-b, + b\ + a.,-a\). (2.10) 

When an operator is applied only to states that are completely in the LLL, it may be represented by a "LLL-projected" 
operator. In the above expression, this means simply dropping the magnetoplasmon operators and a\ that couples 
the LLL states to higher Landau levels. The projected operators then take the form 
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X, -i?,^ + ^(6, + 6j), (2.11) 



y, = i?r + ^(-&. + 6l). (2.12) 
The above operators satisfy the usual LLL commutation relationil 

[Xr,y]]=ll5io- (2.13) 



III. CONDUCTIVITY 



The quantity we will eventually compute is the ac conductivity tensor a^'^{uj), where — x,y. In this section, 
we will show how a^'^{ijj) is derived in terms of the magnetophonon Green's functions. 

The external perturbation due to a spatially uniform electric field with frequency w is H'{t) — eRee^*'^*Eo • u^, 
where = — R^. In the linear response theory, the electric current is given by 

j-(0 = Re^^«(0) (3.1) 

i 

= Re^E / dt'e{t~t'){[ut{tlH\t')]) (3.2) 
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= Re!i^e--X^''(-)K, (3.3) 

where po is the 2D number density of electrons. The brackets (• • •) denote a thermal average for the unperturbed 
Hamiltonian and the zero wave vector response function is defined as 



X^^iu^) = j rft'e^"(*-*'^e(t-i')([<=oW,<=o(i')])- (3.4) 
The linear ac conductivity can be read off from Eq. ( |3.3|) as 



a'"'(c.) = -Im^^X^''(^). (3.5) 

In general, one has to include the polarization field in the calculation and solve for the conductivity self-consistently. 
In a system with uniform background charge, however, the q = component of the polarization is meaningful only 
on the boundaries. In 2D, the boundary effect dies off inversely with the size of the system, and may be ignored in 
the thermodynamic limit. 

The response function may be most easily obtained by computing the Matsubara frequency response function 
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and then analytically continuing 



(3.6) 
(3.7) 



where T r den otes ti me-o rdering in the imaginary time r. Since = r'^ — i?f is a linear combination of bi and b] as 
in Eqs. (2.11) and (2.12), x'^''(iw„) can be written in terms of magnetophonon Green's functions as 
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G-+ G- 

G++ G-i 



The Green's functions are defined as 

^ = - / dre"^-^ ( Tr 
where bb"^ should be interpreted as 6(r)6^(0), etc. 



G-+ G— 
G++ G+- 



1 -i 
1 i 
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(3.8) 



(3.9) 



IV. HAMILTONIAN 



The discussions so far have been quite general and independent of the actual Hamiltonian of the system. In this 
section, we will construct the Hamiltonian in terms of magnetophonon operators. Let us first consider a WC without 
disorder. The kinetic part of the Hamiltonian is given by 



— V + -A 

2m V c / 
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where aj is the magnetoplasmon creation operator in Eq. (|2.7|). The Coulomb interaction takes the form 



He - 



2 K\ri - rj^ 



(4.1) 
(4.2) 

(4.3) 



where k is the dielectric constant. Expanded in terms of the displacement u^, the above equation may be written as 

3" 



where ii,v = x, y, and a is the lattice constant. The dynamic matrix is given by 

^2 / 



ij 



«|R,,|3 



Sfj.1, 



(4.4) 



(4.5) 



where Ry = — Rj. 

Now we take the strong-magnetic-field limit an d pr oject the above Hamiltonian to the LLL. Technically, this means 
rewriting it in terms of ladder operators in Eqs. (2.7) and (2.8), and abandoning any term that is normal-ordered in 

and a, i.e., any term in which all a's are placed to the right of all a^'s. Obviously, the kinetic term is constant after 
the projection, and may be ignored. The Coulomb interaction may be rewritten as 



C2 



E 



4|R,,|3 



b\b, + b,b] + 3 [nlb.b, + nlfblbl 



EE 



2|R,,,|3 



(4.6) 
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where riy = {Rf^ + ii?f^-)/|Ry |. The above expression may also be derived by making use of the analogy between 



the electron system and a spin-lattice system (see Ref. Note that rewriting Eq. (4.4) into Eq. ( |4.6| ), we have 

dropped not only a constant term, but also the last term that is anharmonic. In fact, in order to investigate 
thermal broadening, it is necessary to keep at least one anharmonic term in the Hamiltonian, because a completely 
harmonic Hamiltonian would lead to a theory of free independent eigenmodes and all correlation functions would be 
temperature independent. However, the anharmonic terms in the Coulomb interaction do not contribute to thermal 
shift or broadening of the response function to a spatially uniform external perturbation in a harmonic pinning model. 
This is a highly non-trivial statement, which we discuss in the next section. As a consequence of this, we are led to 
investigate anharmonicity in the other term in the Hamiltonian: the pinning potential. 

For simplicity, we consider only the uniform commensurate pinning model in this work. In reality, of course, 
disorder is more complicated in several ways. First, non-commensurate or non-uniform disorder will break the discrete 
translational symmetry and deform the lattice. Unless the impurity potential is exactly circular at each site, it will 
also mix different angular momentum states. 



However, it was shown in Ref. 18 that at T 0, the uniform commensurate pinning model yields qualitatively the 
same result as more complicated disorder models. Furthermore, experimental data show that the resonance peaks 
keep getting narrower down to T ^ ?ii-<-'p, suggesting that thermal fluctuation affects broadening more than disorder 
does, at least at the observed temperatures. Therefore, we believe this simple disorder model is sufhcient for our 
purpose. 
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FIG. 1. Probability |?/>m(r)p of single particle wave functions in the lowest Landau level. A pit is also schematically drawn 
(a square well at the origin) in order to emphasize that we are considering the case in which the pit size < lo- States with 
m > have little overlap with the pit. 

Fig. shows a schematic picture of the probability function for the three lowest angular momentum states centered 
at an interface pit. Since states with m > have little probability to reside in the pit, they are not affected by it. 
The pinning Hamiltonian is thus appropriately described in terms of magnetophonon operators by 

^P = -^E'5^.^.,o' (4-7) 

i 

where V is the amount of energy an electron pays to move from m ~ state to a higher angular momentum state 
in the absence of Coulomb interaction. This potential is highly anharmonic, and the resulting phonon problem is 
exceedingly difficult to analyze. To overcome this, we adopt an approximate form of the pinning potential. 
The simplest choice that includes anharmonicity is 

Hpi^-jY.^l-blh,)i2~blb,l (4.8) 

i 

which takes the same value as Hp at 6|6i — 0, 1, and 2. At first, the above Hamiltonian looks problematic because 
it has no lower bound. As will be shown below, however, the 1/iV expansion is well defined provided the thermal 
average {b\bi) < 1. It turns out that (fej&i) is indeed small and this approximation gives sensible results in the observed 
temperature range T ^ lOhujp. 

Finally, the total Hamiltonian of our model is given by 

H = Hc2+Hpi. (4.9) 
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A similar Hamiltonian was derived in Ref. |l8| by limiting the Hilbert space to two angular momentum states per 
site and writing the Hamiltonian in a pseudospin representation. Although the higher-angular-momentum states 
should not affect low temperature behavior qualitatively, and the assumption that the many-body wave function 
takes a simple product form breaks down for large m^, the above Hamiltonian has been derived without restricting 
the Hilbert space to a fixed number of states per site. 



V. GENERALIZED KOHN'S THEOREM 



In this section, we will derive a generalized Kohn's theorem, which guarantees that the anharmonicity in the 
Coulomb interaction does not alter the absorption spectrum of a uniform electromagnetic field for a harmonically 
and uniformly pinned model. Similar generalizations have been made in the context of an electron gas in a parabolic 



potential in Refs. |30| and |31|. 

Before we present the proof, let us take a brief look at the original Kohn's theorem. In a magnetic field, the energy 
spectrum of a single free electron is simply given by Landau levels. Now, suppose that there are many electrons 
interacting with one another via Coulomb or any other interaction. Kohn's theorem tells us that the electromagnetic 
absorption spectrum of such a system, as measured through a response to a spatially uniform external field, will be 
exactly the same as that of a single electron. The proof of this theorem is disarmingly simple. Since the electron- 
electron interaction depends only on the relative degrees of freedom, it commutes with the center-of-mass motion. 
Therefore, the observed center-of-mass spectrum is not affected by the interaction at all. 

The big difference in the problem considered in this paper is that there is pinning. However, a similar argument 
can be established in a special case of harmonic uniform pinning model. Suppose that the pinning Hamiltonian takes 
the form 

Hi,,, = VY,blb^ (5.1) 

i 

= Vj2blb^. (5.2) 
q 

Since this Hamiltonian is separable in different momenta q, its eigenstates may be written as a product of the form 

|*) = |*CM>|*rel), (5.3) 

where I^'cm) is an eigenstate of the center-of-mass operator &q=o^q=o and |^rci) of all the other relative degrees of 
freedom. 



Suppose that we calculate the ac conductivity of a non-interacting electron system using the expressions in Sec. HI. 
In order to compute the thermal average, one needs to compute the expectation values over the above eigenstates. 
We make an important observation here that the conductivity depends only on |4'cm), because the operators in Eq. 
(3.4) are functions only of fcq^Q and 6q=o 



Now let us add the Coulomb interaction as in Eq. (O) to the Hamiltonian. As He depends only on relative degrees 



of freedom, it obviously commutes with the center-of-mass displacement 



^E"'- (5-4) 



UCM ^ 

I 

In terms of magnetophonon operators, it implies 

[i?c,foU^q=o]=0. (5.5) 

Therefore, the Coulomb interaction docs not alter the spectrum of I^'cm), although it may (and in general, it does) 
alter that of l^'rci)- Since the ac conductivity is determined by the spectrum of I^'cm) alone, it is unaltered by the 
presence of the Coulomb interaction. This completes the proof of the generalized Kohn's theorem. 

This theorem implies t hat if we use a harmonic pinning potential, the contribution from the anharmonic Coulomb 



interaction terms in Eq. (4.4) to the uniform response function will vanish in all orders of u/a. Therefore, we will 



use a harmonic Hamiltonian for Co ulomb interaction, but an anharmonic one for pinning potential, hence the total 



Hamiltonian is given as in Eq. (4.9). The detailed analysis of the anharmonic pinning potential will be given below, 
using a 1/7V expansion technique. 
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VI. 1/iV EXPANSION 



A expansion technique consists of three steps. First, we make N copies of the magnetophonon operators. 
Second, using a Hubbard-Stratonovich (HS) decoupHng technique, the anharmonic terms are replaced by products 
of a harmonic term and an auxihary HS field. Finally, we expand the magnetophonon Green's function in powers of 
In principle, the expansion is exact if all orders of 1/N are kept. Instead, as is conventionally done, we will stop 
at the first order to get an approximation. Before going into details, it should be noted that N will be set equal to 1 
at the end of the approximation, therefore it is not a perturbative expansion in any sense. However, this technique 
usually does a better job in systematically picking out important processes than many other many-body techniques. 
In fact, the zeroth order approximation, which is a harmonic theory with magnetophonons coupled to the MF value 
of the HS field, already contains the RPA contribution. The first order approximation adds self-energy corrections to 
the Green's functions an d gi ves us a thermal broadening. 

After replication, Eq. ( [4.9D becomes 

N 



1=1 



2N 



(6.1) 



where I, h, and I2 are species indices. Hc2 is the same as in Eq. ( |4.6| ) with different magnetophonon operators 
substituted for each species. There is an extra factor of 1 /N in the last term because it involves two species summations. 
Note that the above equation has been normal-ordered. In this form, we may recast the problem in a path-integral 
form simply by replacing 6 (6^) by a classical field z(t) (z*(t)) in the imaginary time r = it. The Euclidean action is 
given by 



[ li I 



2N 



i L I 



(6.2) 



The first term is the usual Berry phase in the imaginary time. 

In the second step, we introduce a real HS field Qiir) to decouple the quartic term. The procedure is formally 
expressed in a functional integral form as 



exp 



V 
2N 



c / VQi exp 



dT 



V 
2" 



Zli 



where c is a constant. Then the generating functional written in terms of both z and Q takes the form 

Z = j V'^zVQe-^^'', (6.3) 



where 



'S'hs 



dT 



^ z*udrzu + X! Hc2{{zu, 4}) 



(6.4) 



T>^z denotes functional integral over real and imaginary parts of z field. So far the action contains all orders of 1/A^ 
and is therefore exact. In the following subsections, we will Taylor-series expand it in powers of 1/N . 
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A. Mean field theory 



As is the usual practice,Ei we assume that the MF solution of Qi{T) is independent of i and r (Ref. |3^ ). Since the 
crystal momentum q is still a good quantum number, it is useful to Fourier-transform the action. We define 



Zlq 



Zli 



(6.5) 



where g is a shorthand notation for (zct;„,q) and (3 is the inverse temperature. In terms of the above Fourier compo- 
nents, the MF action is given by 



where 



2 



-iq-R,j ^ 



5q 



E 



Rj#0 

The MF Green's functions take the form 

G-+ G—' 
G++ G+- 

1 



(6.6) 

(6.7) 
(6.8) 



w2 + f2 



it^n - /q .9q 

5q it^n - /q 



q Iffql 

We choose the MF value Q in such a way that the MF free energy 

i^MF(Q) - - In / I?2^e-^"''(^) 
has a saddle point. Solving dFup/dQ — 0, we get 



2 p^'^^iF('3)2:*2;q 



-SmfIQ) 



(6.9) 



(6.10) 



(6.11) 



Note that Q above is equal to the thermally averaged value of the number of magnetophonons per site, Tiph = 
^q(6]j6q)/A^s. Q has to be determined self-consistently and is a function of temperature. Note that /q > for all q if 
nph < 1- This implies that the MF theory is stable when there are not too many magnetophonons. The instability in 
the many-magnetophonon regime signals the onset of a "depinning tra nsit ion." In fact, this instability is an artifact 
caused by the fact that our choice of the pinning potential Hp4 in Eq. (4.8) has no lower bound. However, we believe 
there is indeed a depinning transition in real situations, and it is appropriate to regard Upj^ = 1 as the critical number 
of magnetophonons at the transition in this model. As we go beyond the MF theory and include higher orders in 
1/7V, riph decreases due to fluctuations in Qi. 

Since the MF action Smf is quadratic, it may be diagonalized using a Bogoliubov transformation. It is more 
convenient to work in the Hamiltonian representation. The Hamiltonian that reproduces Smf is 



Hmf = /q^l^q - ^ (g^b-qbq + g^b^bl^ 



(6.12) 



We define a new boson annihilation operator 



7q = Uqbq + Wq^Lq- 



(6.13) 
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In order to ensure the boson statistics, it has to satisfy 
By setting 



2 Ve, 



'1 //< 



2 Ve, 



- 1 



we can diagonahze -ffiviF and write it as 



^MF = ^ eq7q7q + COnSt., 



(6.14) 

(6.15) 
(6.16) 

(6.17) 



where = |5q|e"^''. The energy eigenvalues are given by 



= V/q " l5ql' 



(6.18) 



An example of the dispersion relation is shown in Via Jjj , Substituting y = in the above expressions, the dispersion 
relation without pinning, £q oc |qp/^, is recovered. Elll 




\q\a 

FIG. 2. The dispersion relation of a pinned triangular WC is plotted along two directions within the first Brillouin zone as 
specified in the inset. 

Taylor-expanded in small |q|, the parameters and the energy eigenvalue take the form 



/q = (l-g)^ + ^|q|+0[|qn, 



.9q 



2k 

2k Iql 



0[\q\% 



s^^{l~Q)V + —\ci\+0[\ci\% 



which are valid if 



ue 



— |q|«(l-Q)K 



(6.19) 
(6.20) 
(6.21) 

(6.22) 
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One of the authors has obtained a similar dispersion relation at T = 0. The main difference for this finite temperature 
model is that the zeroth order term contains the self-consistent MF value Q. As pointed out in Ref. |l^, the above 
linear dispersion is a direct consequence of the fact that the Coulomb interaction is long-ranged. With this dispersion 
relation, the magnetophonon density of states vanishes at the lower band edge, and thermal broadening will be 
substantially suppressed as will be shown below. This is not the case for a short-range interaction. For example, a 
screened interaction would have yielded a dispersion that is quadratic at small wave vectors. In that case, the density 
of states jumps to a finite value at the band edge. In a separate test, we have found that thermal broadening is indeed 
several orders of magnitude greater for a screened interaction than for an unscreened one. 

Since /q and Eq are functions of Q, the saddle point solution for Q has to be determined self-consistently. More 
specifically, it has to satisfy 



]^ E I"ql'(7/q7/q> + bql'(7i,-q7/,_q) 

* q 



A /q 1 



2 7 £q 2 



(6.23) 
(6.24) 
(6.25) 



where nB{s) = XjieP^ — 1) is the Bose function. 

In the MF theory, the uniform field absorption spectrum will have a delta function peak at TiWp 



eq=o- Since 



Eq^o = V^(l ~ Q) and Q increases with temperature, the thermal reduction of the peak frequency observed in 
experiments is explained qualitatively already at the MF level. Physical interpretation is also simple. At high 
temperatures, more magnetophonons are thermally created. This implies that the average number of electrons that 
stay in the pits decreases, leading to a drop in the spatially averaged pinning potential. 

In order to explain thermal broadening, however, one has to go beyond the MF theory. We will compute the first 
order corrections in the expansion in the next subsection. 



B. Order IjN corrections 



In this subsection, we will consider fiuctuations of the HS field about its MF solution. The fiuctuation fields are 
defined at each site as 



and its Fourier-transformed field is 



Since riir) is real, r_g — r*. In terms of r^, Eq. ( |6.4D may be rewritten as 



(6.26) 
(6.27) 



= -f E E [ ^'-9 ] (^MF + U),,' 
; qq' 

(3NV 



Zlq' 

<-q' 



where 



i^n - /q 



.9q 


-iuJn - /q _ 


1 


■ 




Vrq^g 



Sqqi , 



(6.28) 

(6.29) 
(6.30) 



Gmf is the MF magnetophonon Green's function matrix and U describes the coupling between a magnetophonon and 
the HS field r, or a vertex in the Feynman diagrams. The "bare" HS propagator can be read off from 5hs a-nd is 
given by 
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NV 



Fig. |^(a) shows the building blocks of the Feynman diagrams that will be frequently used below. 

- C 



> 








—>- 







''MF 



~ "-^MF 

= Do 



— 
V 




FIG. 3. (a) The Feynman diagrams are drawn for MF magnetophonon Green's functions Gmf (straight line), Hub- 
bard-Stratonovich propagator Do (curly line), and the vertex V/y/W^. The arrows on both sides of a vertex must be the 
same. Note that there are anomalous Green's functions that have two arrows pointing opposite directions, because the Hamil- 
tonian does not conserve the magnet ophon on number, (b) Examples of loops with vertices. Each loop corresponds to a 
summand of the n-summation in Eq. (3.34). n — 2,3, and 4 respectively. A straight line without arrows denote a sum of all 
arrowed lines that are permitted by the "directi on-co nservation" rule at each vertex. The first diagram represents the RPA 
propagator. Since there is no n = 1 term in Eq. (6.34), the Feynman rules do not allow a loop to have only one vertex. 



Note that Dq is proportional to 1/iV, so each curly line raises the order of an overall diagram by one in the 1/iV 
expansion. On the other hand, a magnetophonon loop [see Fig. ||(b)] contributes order TV due to a summation over 
an internal species-index. In order to find corrections to the HS propagator, it is convenient to integrate out the 
magnetophonon degree of freedom. The effective action for the HS field is given by 



S,s[r] = - In jv 



N 



^--tr In [/3(G 



N 



-1 

MF 



U)] 




(6.32) 



(6.33) 



tr In 



f3NV 



El 

9 



E: 

n=2 



-tr(GMFU)" 



(6.34) 



Going from Eq. ( 6.33| ) to Eq. ( 6.34 ), we have Taylor-series expanded the logari thmi c function about /3G^p. Unim- 
portant constants were omitted at ever y step. Note that the last term in Eq. ( 6.33| ) exactly cancels out the n = 1 
term that would have been in Eq. ( ^.34| ), because Q is chosen at a saddle point. Gmf and U are considered as entities 
in a product space of a 2 x 2 matrix and a matrix labeled by frequency and wave- number indices qq' . Thus, a trace 
operator sums not only over the 2x2 matrix indices, but also over q. 

Let us take a closer look at the last expression ( ^.34 ). The first term is the MF free energy and the last term is 
the action of free independent HS fields r^. Each of the remaining terms is represented by a loop with n vertices as 
in Fig. |(b). 



FIG. 4. Dyson equation for the RPA Hubbard-Stratonovich propagator, 
for the RPA propagator. 




Thick wiggly lines denote self-consistent solution 
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Let us now search for the fowest order correction to the HS propagator . The re is no correction of the zeroth order 
or lower, and the only 1/N correction conies from the n = 2 term of Eq. ( 6.34 ). As can be seen in the first Feynman 
diagram in Fig. ||(b), this term is nothing but the RPA correction to the HS propagator. We explicitly write it as 



— trGMpUGMpU = — - — QMFiqjp 



ql 1 



(6.35) 



where £>mf('z) is the MF magnetophonon density-density correlation function defined as 



£'mf('7) 



2(3K 



j!— ^trGMF(9' - q)GMF{q') 



MF 



(6.36) 



(6.37) 



where (• • •)mf denotes an average in the MF theory. Detailed calculations of (?mf('z) will be given in Appendix A. 
The RPA "dressed" propagator is defined through the Dyson equation 



1 



D^' +NV^gMF{q) 
1 

~ NV{l-VQMF{q))' 



(6.38) 
(6.39) 



Since both the bare propagator Do and the RPA correction are of order the above propagator is also of order 

A thick curly line will be used for the Feynman diagram of the dressed propagator as in Fig. 0. 



op op „ 



q-q 



op 



+ 




FIG. 5. Two self-energy diagrams of order 1/N. The indices ct and p take either -I- or — . They also determine the direction 
of arrows in the diagram of Ea. A — (-I-) sign means annihilation (creation) of a magnetophonon, so the corresponding arrow 
points towards (away from) the nearby vertex. 

We now have all the ingredients to calculate the correction to the magnetophonon Green's function. Keeping 
in mind that each HS propagator is of order 1/N and each magnetophonon loop is of order N, it is straightforward to 
count the order of any Feynman diagram. There are two self-energy diagrams of order 1/N as shown in Fig. ^ They 
are written in terms of MF Green's functions and the HS propagator D as 



K'iq) 



NV^ 



1' 

D{Q)Y,D{-q'-q") 

q'q" 



(6.40) 



>^<G^W)[{GM%{q)f + GMF{q)Gl,^{q) 



Gl,^{q")Gt^^{q')G^l{q') + C.C 



(6.41) 



These expressions can be partially evaluated analytically via frequency summations .c3 More details may be found in 
Appendix B. 

The self-energy correction is incorporated into the Green's functions via the Dyson equation 



(6.42) 
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where 



1 = 


'10" 




1 



(6.43) 



Combining Eq. jel^ ) with Eqs. (fj), (js^), and (U), we can finaUy compute the ac conductivity to the order 1/N. 
Using the final expressions in Appendix B, we have numerically computed the self-energy diagrams and obtained the 
conductivity as a function of T. The result is presented in the next section. 



VII. RESULT 

Before we present the result, it is useful to review some important energy scales. For the purpose of comparison, 
we will use the data in Ref. ^. We find 

— — ^ 16 meV, 
kIo 



TlLO, 



na 
huj„ 



> 3 meV, 
~ 2 meV, 

< 60 ^eV, 
- 5 /leV. 



(7.1) 



The temperature range for which the line width was clearly discernible was 3 /ieV < ksT < 18 ^eV. 
For efficiency in numerical computations, we have made the following isotropic approximations 



/q = /(|q|), 
5q = .9(|q|) 
£q = e(|q|)- 



iqy 



(7.2) 



We did, however, keep the shape of the first Brillouin zone to be hexagonal. According to our numerical tests, these 
approximations break down as the Brillouin zone boundaries are approached, i.e., when |q| ~ 7r/a, but they are very 
accurate for |q| < 1/a. For example, when V ~ fiLOp, the energy Cq is already an order of magnitude greater than hup 
at |q| — 1/a (see Fig. ||), but the error Ae is still less than 4%. Therefore, this approximation should be quantitatively 
reliable if T is not too much greater than hujp. 



a 
a 



^ I — ' — I — ' — r 

8 -aotonooo 

-r_i_ w I 

7 - 

6 - 

5 - 

4 - 

3 - 

2 - 

1 - 



o expansion 
+ Experiment 



FIG. 6. Peak frequency as a function of temperature, 
(Ref. H), but there is a good qualitative agreement. 



100 200 300 400 500 600 
T (mK) 

Our results are overestimated compared to the experimental data 
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In order to compare our result with the experimental data, we have chosen the bare pinning potential V in such a 
way that the resulting peak frequency tUp agrees with the experiment at the lowest measured temperature. The peak 
position data are shown in Fig. ^, along with experimental data from Ref. ||. Although our results are systematically 
overestimated, the qualitative agreement is rather good. 



0.01 



0.0001 




'^"^^0 100 1000 

T (mK) 

FIG. 7. The half-width-half-maximum line width is plotted against temperature in log-log scale. The peaks appear to be 
extremely sharp especially in the low temperature regime. However, the width increases very rapidly with temperature. 

The results for line width Auj are surprisingly small (Fig. |^) . Compared to those from experimental data, our results 
are several orders of magnitude smaller at the lowest T. This is especially surprising because T ~ tiujp. Naively, one 
would expect that the line width should be of the same order of magnitude as T, due to thermal broadening. Although 
much greater than our results, experiments also confirm that this is clearly not the case. In Refs. |^ and ||, Alu is 
almost one order of magnitude smaller than ojp. In Ref. ^ the quality factor Q = ujp/Aoj even exceeds 30. There are 
also other indications that the peaks may be actually even sharper.! 



The unusually sharp peak may be rather well understood within our model. As discussed in Sec. VIA, due to the 
long-range Coulomb interaction, the dispersion starts out linearly in |q|. The density of states thus increases linearly 
from zero. Explicitly, the density of states per unit cell takes the form 



AA(e) 




[e - (1 - Q)V] if e > (1 - Q)V, 
if £ < (1 - Q)V. 



(7.3) 



A rough estimate of the total number of "accessible states" per site at temperature T would be 



Oc / de M{s), 



(7.4) 



where ac is the unit cell area. Using the parameters in Ref. ||, we find that there are merely 2 x 10~^ accessible 
states per site at T = 180 mK, which is the highest temperature for which the line width was measured. The physical 
interpretation of this may be given as the following. At least in the semiclassical level, our model is essentially that of 
a crystal in which each electron is attached to a lattice site by a short-range binding potential. The overall uniform 
motion (q = 0) determines the peak frequency. Broadening comes about from thermal and quantum fluctuations 
deforming the lattice and making the electrons move relative to one another. However, the above analysis tells us 
that a WC in a strong magnetic field is so rigid that it is hardly deformed even at temperatures greater than the 
pinning potential. 

In reality, the external ac field has a finite wave length of order ^ 30 /im ~ lO^a (Ref. 0). Accordingly, we have 
performed a similar analysis for the response at the finite wave vector, but we found no qualitative change. 

Defects in a WC such as interstitials, vacancies, dislocations, etc. may also affect broadening. It is certainly true 
that their presence will soften the lattice and subsequently enhance magnetophonon excitations. However, unless one 
is too close to the melting temperature of the WC, which is not the case in the experiments, their effects may be 
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mostly taken into account through renormahzation of the lattice stiffness, or dynamic matrix in Eq. (4.5). This should 
not affect our results qualitatively. 

There may be other sources of broadening, too. Other low energy modes, which are-not included in this work, 
may be involved in the broadening. One possible candidate is the edge states of the WC,L3 the analysis of which will 
be given elsewhere. Another possible source of broadening is extrinsic low energy modes outside the WC. When an 
experiment is performed, there are many external degrees of freedom that may be coupled to the sample. If some 
of them are hard to remove and thus are left undetected to influence the data, the observed peak will be certainly 
broader than theoretically predicted. l— . 

It has also been suggested that higher Landau level mixing might be important .E3 It is certainly true that there 
can be second or higher order scattering processes that use higher Landau levels as virtual states. However, these 
processes should be suppressed as a power {up/ucY with a > 1. Since Wp/wc 10"'^, Aw as a result of such processes 
would be still too small to explain currently available experimental data. 



VIII. CONCLUSION 



In this work, we studied the thermal broadening of the electromagnetic absorption resonance of a magnetically 
induced WC. In the strong magnetic field limit, the low-energy collective modes that are coupled to the spatially uni- 
form ac electric field are magnetophonons. Assuming that electrons are closely bound to lattice sites, magnetophonon 
creation and annihilation operators may be constructed out of displacements = r; — R^. 

We showed that the Hamiltonian could be divided into harmonic and anharmonic parts in terms of magnetophonon 
operators. The harmonic part, when taken alone, describes independent magnetophonons and produces a delta 
function peak in the spectrum. The anharmonic part introduces magnetophonon interactions. It will not only 
renormalize the peak position, but also broaden it by mixing different magnetophonon modes. Anharmonicity in our 
model comes from two sources: the Coulomb interaction and the pinning potential. Since the Coulomb interaction 
depends only on the relative coordinates, it is completely decoupled from the center-of-mass degrees of freedom. This 
leads to the derivation of a generalized Kohn's theorem that asserts that the Coulomb interaction cannot shift or 
broaden an ac conductivity resonance peak even in the presence of a uniform harmonic pinning potential. 

Analysis of the magnetophonon interactions in the pinning potential was performed using a 1/N expansion tech- 
nique. This technique provides a systematic way of summing up important diagrams and it captures important pieces 
of physics in the low order solutions. For example, we found that the zeroth order MF solution could account for the 
decreasing peak frequency LOp as a function of temperature. 

Thermal broadening appears in the 1/A^ self-energy corrections. However, the line width Aw is found to be many 
orders of magnitude smaller than ujp and T. The reason for this lies in the magnetophonon dispersion that is linear 
and steep in the low energy limit. This translates into a small density of states for collective excitations. Consequently 
the broadening is substantially suppressed. „ „ 

Our result of peak frequency agrees qualitatively with recent experiments. Q'u The line width result is also qualita- 
tively consistent in that it is much less than ujp and T. Quantitatively, our result of Aw appears to be much smaller 
than published results, although there is continuing experimental work on sorting out just how narrow the intrinsic 
line width really is.o Other possible sources of broadening, such as edge states and extrinsic low energy modes, may 
also be responsible for observed line width. 
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APPENDIX A: DENSITY-DENSITY CORRELATION QMF(q) 

In this Appendix, we will present detailed calculations of the density-density correlation function qmy{<1)- Eq- ( 6.36| ) 
can be written as 
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9mf 



f3N, 



(Al) 



where we have used 

Gl^A-Q) - G£[p(g) (A2) 
to get the second line. Diagrammatically, the last two terms are represented by two arrowed diagrams as in Fig. ||. 






FIG. 8. Two arrowed diagrams for the density-density correlation function £iMF(fi'). 
Using the Green's functions in Eq. ( |6.9| ), we may rewrite the above equation more explicitly as 



QMFiq) 



—y 

q 



{n, 



(eq'-q - £ q')(l +^q'~q.q') 
(Sq'-q ~ £q')^ 



q'-q "q'; /. 



+ (n, 



q'-q 



1)- 



(£q'-q + £q')(l-6 



(iw„)2 - (eq'_q + Eq'Y 



(A3) 
, (A4) 



where riq' is a shorthaad notation of the Bose function nsisq'). In the second line, we have used the usual frequency 
summation technique.EJ 

g{q) has many simple poles in the complex to plane on the real axis. In the thermodynamic limit, they become 
branch cuts. Therefore, summations involving g{q) must be performed with care, as will be demonstrated in the 
calculations of self energy corrections in the next Appendix. Note that g{q) is analytic and real at w = 0, so that 
the branch cuts are divided into two: one on the positive and the other on the negative axis. Both branch cuts are 
bounded because Sq'-q ± Sq' is finite. 

Eventually, we will analytically continue itOn lo + iO^ to get a retarded correlation function. Since g{~-q) — g{q) 
and g{q) is analytic except on the branch cuts, it has the following properties: 



Re g{—uj + iO^, q) = Re g{iLj + iO^ , q), 
Im g{'-LLj + iO^, q) ~ — Im g{LLi + iO^ , q). 



(A5) 



which will become useful in the next Appendix. Finally, g{uj + iO+,q) was numerically computed using an isotropic 
approximation as explained in Sec. VII. 



APPENDIX B: 1/N SELF ENERGY CORRECTIONS AND Es 

In this Appendix, we will calculate the self energy corrections to the magnetophonon Green's functions as in 
Fig. ^. Since the bare HS propagator Dq is a constant, it is convenient to separate it from the rest of the RPA dressed 
propagator D. The remaining part 



SD{q) = D{q) - Do 



N[Vgo{q) - 1] 



(Bl) 
(B2) 



vanishes as w ^ in the large lu limit, so all following frequency sums that contain (5-D(q) converge without introducing 
cumbersome convergence factors such as e^*'^'^ . 



The first self-energy correction is given in Eq. ( 6.4C ). The contribution of the bare propagator Dq is simply 



(B3) 
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The summation of MF Green's function is defined with an appropriate convergence factor that ensures normal- 
ordering. It is then straightforward to show 



EGmf(9')--/3E(VV)^ 

q' q' 



(B4) 
(B5) 
(B6) 
(B7) 



q' 



where rian is the anomalous magnetophonon number per site. Due to the six-fold symmetry of the lattice, however, 
nan = 0. We therefore get 



(0 



— h 

aO 




aO 


- y++ 


X 





N ' 



(B8) 
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> 
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FIG. 9. Contour is deformed in the complex u). Simple poles are denoted by x and branch cuts by zigzag lines. The 
separation between two branch cuts is exaggerated in order to show that the pole at the origin is treated separately. In fact, 
the branch cuts asymptotically touch the origin. 

Now let us turn our attention to the remaining part. It needs to be computed with more care, since SD{q') has 
branch cuts on t he real u axis, that arise from the density-density correlation function g{q'). Using the Green's 
functions in Eq. (3.£), we may write 



N(3N, 



NfiN, 



(B9) 
(BIO) 



where Mq' and v^^i are defined in Eqs. ( |6.15| ) and (3T£). As in a standard frequency summation technique J^l we may 
replace the sum over iuo'^ by a contour integral on the complex lo' plane with simple poles on the imaginary axis. 
After the contour is deformed as in Fig. ^, the new expression takes the form 



^ [5D{iuJn - £q',q- q')|uq'P + 5D{iuJn + eq',q- qOlvH "s(eq') 



ttNN. 



E 

q' 



dLj'[Im6D{u;' + iO+ q')] 



|Uq'P"i3(w') IWq'pnB(-w') 



icUn + Uj' — Ec, 



lUJn + IjJ' + £c 



(Bll) 



where (jj' = is excluded from the domain of the integral 

, -r 



0- . 



(B12) 
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Although UBiio') diverges as w' 0, the above integral is well defined, because lmSD{uj' + iQ~^, q — q') vanishes as 
a;' = is approached. Finally, analytic continuation is performed to give 

sj:-+{lu + ?:o+, q) = ^Y. [^^(^ - + q - q')l"qf + SOicj + e^, + iO+, q - qOKv P] neie^.) 



q' 



ttNN. 



LU + LU — Sq' LU + LU' + Sq' 



J2^J dcj'ImSD{Lj' + iO+,q - q') 



'^^inSD{Lo - Sq' +iO+,q- q')|uq'pn_B(eq' ~ lu) 
q' 

+lin5D{uj + £q' + iO+, q - q')\vq'fnB{eq' + uj) 



(B13) 



where Pj denotes a Cauchy principal integral. Eq. (A5) was used in the last line. The above quantity may be 
computed using the numerical solution of p(a; + iO^, q). The other JSq's can be calculated in a similar way. The final 
form of is the same the above equation, except that |tiq'P and |fq'P are switched. ST,~~ is obtained 

by replacing both Iwq'P and |fq'P by u*,Vqi (uqtv*,). In the isotropic approximation, the angular sum of 



U„iV, 



q'Vq' 



2e, 



q' 



« - i<lv 



vanishes, so 



-+ _ 



0. 



(B14) 



(B15) 



The other 1/iV correction, E^, may be calculated in a similar manner. First, we compute the contribution from the 
bare HS propagator. Substituting Dq for D{—q' — q"), Eq. (6.41) becomes 



-'bO 



13^ 



q'q" ^ 



Using Eqs. (B4) through (B7) along with 



(B16) 



we get 



E(Gmf('z'))'=E 



|Uq'| 



ILOn - En' lUJn + £. 



q' 



q 

4\ 



{\uq,\' + \v^t)n'B{eq,)-2\uq^^\v^\- 



2e, 



q' 



lUJn - £q' l^^n + £q' 



-/^EKI'IvP 

q' 



n'g{eq,) 



2na(gq>) + l 



2e 



q' 



[Vq{Q) - l]NN, 



E 



ffq'P 2 nB(£q') + 1 /q' . 

^2^. ^"s(eq') 



q' 



(B17) 



(B18) 



(B19) 



The derivative of the Bose function satisfies 7i^(£q') = —l3nB{£q')[nB{£q') + 1]- 
For the remaining part, we substitute SD{—q' — q") for D{—q' — q"). We get 



/VT/4 

q'q" 



SD{-q' - q") 



\u'u"* + v'v"*\^ + Iff) \u'v" + v'u"\^ + \v'\^) 

+ £')' ^ - 

2(u'u"* + t>'w"*)(u'i;" + + c.c. 



(B20) 
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Shorthand notations are defined as u' = u^', u" = u^", etc. In order to avoid the branch cuts in SD, we make a 
change of variables icu!^ — > — — ii^n and perform the frequency sum over zw". The result takes the form 



HBie") - nsis') n'sie') 



|uV' + fV'|2f|u'|2 + |f'P 



{iu'^ - e' + e")^ iuj'„-e' + e" 
neie") + neie') + 1 n's(e') 



{itu'^ + e' + e")^ iuj'^ + e' + e" 
{u'u"* + v'v"*){u'v" + v'u"y*v'* + c.c. \nB{e") - UBie') UBie") + ns(e') + 1 



iujL -e' + e" 



iujL + £' + e" 



(B21) 



The last frequency sum X^ic^' performed using a similar contour deformation technique as in Fig. |^. Using 

, 5D{uj' + iQ+)nB{uj') 



^^(^-"1 ^ ^Im / d.' 



lUJ' — e TT 



= -Pi djjj' 

V /^(-"L = ^Im / d.' 



w' — e 



1 {i^'n - 



— 5D{uj' + iQ+)nB{uj') 
oui' 



/?Re(5D(e + iO+)ns(£), 
1 



(B22) 

(B23) 
(B24) 



Eq. ( B21 ) becomes 



[Vg{0) - 1]N! 



C2 



■^,5D{J,-c{-c(')nB{Lo') 



C3 



+ iO+ uo' + e' + e" + iO+ 

C4 



w' - e' + e" + iO+ w' + e' + e" + iO+ 



(B25) 



where 



I / //* / //*l2 /I n2 I n2\ / / /N iu'u"* +V'v"*){u'v" +V'u")u'*v'* + C.C. 



[uBie') ~ nB{e")], (B26) 



, , , ,,,2. 1 „2 , I n2N , . n , (uV* +i;V'*)(uV' + t;'u")u'*'y'* +C.C. 
C2 = -|wV + -y'-u"!^ [\u'\^ + Iw'l^) n^(e') + -^^ ; 



[nB(e')+«s (£") + !], (B27) 

(B28) 



C4 = \u'v" + fVf (|uf + ) [ns(£') + nB(e") + 1]. (B29) 
Finally, combining Eqs. (p8|), ( |B13 ), ( B1£ ), and (B25), the order self-energy correction is given by 

S7 - + <5S7 + (Sbo + 5^b)5.p. (B30) 
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